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I.  INTRODUCTION 

A.         BACKGROUND 

Propagation  of  radiowaves  over  rough  ocean  surfaces  and  in  an  inhomogeneous 
atmosphere  is  a  topic  of  particular  interest  to  the  Navy.  The  effects  of  this  type  of 
environment  on  the  strength  of  a  signal  at  a  given  distance  downrange  from  a  transmitter  are 
significant,  and  play  a  critical  role  in  determining  if  a  communication  link  can  be 
successfully  established,  or  a  target  can  be  detected  by  radar.  This  type  of  radiowave 
propagation  is  governed  by  a  modified  Helmholtz  equation  which  is  an  elliptic  partial 
differential  equation.  Solving  this  type  of  equation  is  computationally  very  intensive  and 
becomes  impractical  when  addressing  typical  propagation  problems.  The  Helmholtz 
equation  can  be  approximated  by  a  parabolic  equation  (PE)  if  only  one  way  propagation  is 
considered,  which  is  a  reasonable  assumption  for  most  typical  propagation  problems.  The 
advantage  of  using  the  parabolic  equation  approach  is  that  a  solution  can  be  obtained  much 
more  efficiently  through  the  use  of  a  range  stepping  technique.  Kuttler  and  Dockery  [3] 
discuss  the  basic  idea  and  various  approximations  involved  in  the  development  of  the  PE. 
This  technique  makes  it  possible  to  easily  estimate  propagation  losses  several  hundreds  of 
kilometers  downrange  with  antenna  heights  up  to  several  hundreds  of  meters  for 
frequencies  through  the  Super  High  Frequency  (SHF)  band.  By  only  considering  waves 
propagating  in  the  forward  direction,  the  strength  of  a  signal  at  a  given  location  downrange 
will  be  determined  by  direct  and  reflected  (from  the  ocean  surface)  waves.  Interference 
between  these  two  waves  will  result  in  reflection  multipath  fading  which  can  result  in  high 
gain  or  severe  loss  of  signal.  The  effects  of  reflection  multipath  fading  will  depend  greatly 
upon  the  surface  height  deviations  of  the  ocean  surface  which  will  be  determined  by  wind 
speed. 


B.         OBJECTIVE 

In  this  thesis  the  split-step  PE  algorithm  is  used  to  predict  propagation  of  vertically 
polarized  radiowaves  over  rough  ocean  surfaces  in  the  presence  of  an  inhomogeneous 
atmosphere.  This  thesis  builds  upon  work  previously  done  by  Janaswamy  [1]  in  this  area 
which  addressed  horizontal  polarization.  In  this  approach  the  Fourier  transform  pair  used  in 
the  split-step  PE  algorithm  is  modified  to  accommodate  the  effects  of  sea  surface  roughness 
directly  into  the  formulation.  The  idea  behind  this  approach  will  be  to  cast  the  transform 
equations  in  terms  of  incident  and  reflected  waves,  and  then  use  the  rough  surface  reduction 
factor  available  for  plane  waves,  according  to  Miller  [2],  directly  in  the  spectral  domain. 
Chapter  II  presents  the  derivation  and  formulation  of  the  modified  Fourier  transform  pair 
for  the  split-step  PE  algorithm.  Chapter  III  details  the  generation  of  the  numerical  procedure 
for  solving  the  parabolic  equation.  In  Chapter  IV  the  performance  of  the  numerical  solution 
is  examined  by  varying  several  important  parameters  and  observing  the  effects  on 
radiowave  propagation.  The  effects  of  changing  parameters  such  as  step  size,  atmospheric 
refractivity  profile,  number  of  points  for  Fast  Fourier  Transform  (FFTs),  frequency,  and 
wind  speed  will  be  investigated.  Recommendations  and  conclusions  are  presented  in 
Chapter  V. 


II.  FORMULATION 

In  this  chapter  we  present  the  theory  governing  the  parabolic  equation  and  the 
derivation  of  the  Fourier  transform  pair  to  be  used  with  the  split-step  algorithm.  This 
derivation  parallels  that  given  by  Janaswamy  [1]  which  considered  the  same  propagation 
problem  discussed  here,  but  for  horizontally  polarized  radiowaves.  This  chapter  reiterates 
many  of  the  formulas  from  that  report  but  modifies  them  as  appropriate  to  derive  correct 
formulas  for  the  case  of  vertically  polarized  waves. 

Figure  1  illustrates  the  basic  geometry  of  the  problem  we  are  investigating.  We 
consider  the  source  at  an  initial  range  (x  =  0)  and  height  (z  =  z,)  to  be  an  omnidirectional 
point  source.  Given  the  position  of  the  receiver,  frequency  of  operation,  refractive  index 
profile  of  the  atmosphere,  wind  speed,  and  the  ground  constants  (sr  and  a),  we  wish  to 
determine  the  signal  strength  along  the  path  between  the  transmitter  and  receiver. 
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Figure  1.  Source  Producing  Fields  Over  Rough  Ocean  Surface 


To  solve  this  propagation  problem  the  standard  parabolic  equation  as  defined  below 
by  Kuttler  and  Dockery  [3]  is  used: 


d2u(x,z)  du(x,z)  2r    /       \       I   /       \ 

— —^ —  +  2ik0  — — —  +  2k0  \m[x,z)  -  \\u[x,z)  =  0 


(1) 


/      x       Irsin0      ,     v 
where    u{x,z)=   l—j — -rH^{r,v)    refers   to    a   spherical    coordinate   reference    system. 

H^(r,Q)  is  the  (j) -component  of  the  magnetic  field.  Equation  (1)  assumes  an  e_"°'time 
dependence.  The  coordinate  system  used  in  this  thesis  is  shown  in  Figure  2  below. 


►  Y 


Figure  2.  Earth  Centered  Spherical  Geometry 


In  this  coordinate  system  the  transmitter  is  located  at  0  =  0  and  r  =  ae  +  zt  where  ae  is  the 
radius  of  the  earth  and  zt  is  the  height  of  the  transmitter  above  the  surface  of  the  earth.  The 
free  space  wave  number  is  given  by  k0  =  co^/e^"  ,x  is  the  range  axis,  z  is  the  height  axis 
and  n\x,z)  is  the  modified  refractive  index  under  earth  flattened  conditions  as  given  by 

Kuttler  and  Dockery  [3]  and  is  equal  to:  m(x,z)  = ,  with  n  being  the  actual  refractive 

index. 

By  examining  (1)  we  can  see  that  the  highest  derivative  in  terms  of  the  range,  x,  is 
of  the  first  order  for  the  parabolic  case,  unlike  a  Helmholtz  equation  which  would  be  of 
second  order.  This  desirable  characteristic  of  the  parabolic  equation  allows  us  to  use  (1)  to 
solve  for  a  field  at  a  given  range  based  upon  a  known  field  at  a  previous  range.  If  the  range 
step  size,  or  Ax,  is  kept  reasonably  small  and  the  refractive  index  varies  slowly,  the  field  at  a 
new  range  can  be  computed  from  the  field  at  a  previous  range  by  using  the  split-step 
algorithm  by  Tappert  [4]  below: 


u{x,z)  =  e'^m-1)Ax¥-l< 


-ip26x 


e2k 


°  f[w(x0,z)] 


(2) 


Later  in  this  thesis  we  will  modify  this  equation  to  accommodate  refractive  indices  which 
have  steep  gradients.  F  represents  the  forward  Fourier  transform  operator  corresponding  to 
an  appropriate  spectral  decomposition  of  the  field  in  the  vertical  direction,  and  F"1  is  its 
inverse. 

To  develop  the  form  for  the  Fourier  operators  to  be  used  in  the  split-step  algorithm 
for  rough  ocean  surfaces  it  is  logical  to  first  solve  for  the  smooth  surface  case,  and  then 
modify  this  result  so  it  can  be  applied  to  a  rough  ocean  surface.  In  our  formulation  the  ocean 
surface  shows  up  as  a  boundary  condition  for  which  the  parabolic  equation  must  be  solved. 
For  the  smooth  ocean  surface  case  the  parabolic  equation  as  given  at  (1)  must  be  solved 
subject  to  the  impedance  boundary: 


du(x,0) 


dz 


+ 


a0w(x,0)  =  0  ,  (3) 


where  a0  = - .  The  complex  dielectric  constant,  src  is  defined  as: 


re 


Zrc=Zr+i°r,  (4) 

where,  ar  is  the  relative  conductivity  of  the  earth  given  by: 

<5r  = ,    co  =  angular  frequency  (rads  /  s)  (5) 

cos0 

The  boundary  condition  given  at  (3)  is  valid  only  for  |src|  »  1 ,  which  for  ocean  surfaces  is 
easily  met.  To  solve  the  standard  parabolic  equation  given  at  (1)  subject  to  the  boundary 
condition  at  (4),  the  following  mixed  Fourier  transform  pair  by  Kuttler  and  Dockery  is  used 
[3]: 

oo 

u(x,p)  =  ¥s(u)  =  jw(x,z)[a0  sm(pz)  -  pcos(pz)jdz  ,  (6a) 


(      \     v-U~\     2TV       \^0sm(pz)-pcos(pz) 

u[x,z)  =  ¥iu)  =  —  \u[x,p) 2 1 ap  +  S\x)e    * 

« o  ao+P 


(6b) 


The  subscript,  s,  stands  for  the  smooth  surface  case  and  the  last  term  in  (6b),  which 
represents  a  surface  wave,  decays  with  range  and  height  and  can  be  ignored  for  frequencies 
over  10  MHz  according  to  Kuttler  and  Dockery  [3],  which  will  be  the  case  for  our  work 
here.  To  better  illustrate  the  fact  this  formula  represents  plane  waves  which  are  traveling 
towards  and  away  from  the  boundary,  it  is  rewritten  in  the  mathematically  equivalent  form 
shown  below: 


1 


u 


(x,p)  =  ¥s(u)  =  \u(x,z)eipzdz  +  —r^-  ju(x,z)e'ipzdz  ,  (7a) 

0  *  s[P)  0 


2  °° 
u(x,z)  =  F,(k)  =  -  \u{x,p)[e"pz  +  r,(p)e*»]#  ,  (7b) 

/  x      p  +  ia0 

where, Ts{p)  = ; —  ,p  >  0,  represents  the  plane  wave  reflection  coefficient  for  smooth 

p  —  i(x0 

earth.  From  examination  of  (7b)  it  is  apparent  that  the  field  is  comprised  of  a  wave  traveling 
towards  the  boundary  at  z  =  0,  given  by  u(x, p)e~,pz ,  and  one  away  from  the  boundary  at  z  = 

0,  given  by  Ys\p)u\x,p)eipz .  It  is  evident  that  the  reflected  wave  is  equal  to  the  incident 
wave  multiplied  by  the  reflection  coefficient,  which  makes  intuitive  sense. 

Now  that  Fourier  operators  have  been  defined  for  the  split-step  algorithm  in  the  case 
of  a  smooth  ocean  surface,  we  want  to  extend  our  derivation  to  include  a  boundary  which 
consists  of  a  rough  ocean  surface.  Previous  investigation  of  the  problem  of  plane  wave 
reflection  from  rough  surfaces  has  resulted  in  the  Kirchoff  approximation  and  the  concept 
of  a  rough  surface  reduction  factor  as  given  by  Beckmann  and  Spizzichino  [5].  To 
accommodate  surface  roughness  for  our  problem  we  simply  modify  the  smooth  surface 
reflection  coefficient,  Ts(p),  by  multiplying  it  with  a  rough  surface  reduction  factor  (defined 
below)  to  yield  a  rough  surface  reflection  coefficient,  TT(p). 

rr{p)  =  Po{p,vhYs{p)  ■  (8) 

The  term,  p0(p;aj,),  is  the  rough  surface  reduction  factor  and  has  been  derived  by  Miller  [2] 
as: 

p0(p,Gh)  =  e-2p2°h0(2pWh).  (9) 

The  I0(-)  term  in  (9)  is  the  modified  Bessel  function  of  the  first  kind  of  order  zero  and  ah 
represents  the  rms  wave  height  deviation  as  determined  by  wind  speed  according  to  the 
formula: 


ah  =  0.005  l^i2  ,  (10) 

where    \i  =  wind  speed  (m/s).  A  simpler  form  for  the  rough  surface  reduction  factor  has 
been  recommended  in  CCIR  report  1008-1  [6]: 

J 

^32X -2  +  yl(32X)2  -IX +  9 


Po{P'°h)=     I  r =  (11) 


where  X  =  2p2a2h . 

We  have  now  presented  the  necessary  material  to  allow  us  to  modify  the  Fourier 
operators  such  that  the  split-step  algorithm  can  be  used  to  predict  propagation  of  vertically 
polarized  waves  over  rough  oceans  surfaces.  The  field  equation  given  for  smooth  surfaces  at 
(7b)  now  becomes  as  follows  for  rough  surfaces: 

1    °° 
w(x,z)  =  F;1^*,/?)]  =  —  \u{x,p\e-,pz  +  Tre,pz]dp  ,  (12) 

where,  Fr_1 ,  represents  the  inverse  Fourier  operator  for  rough  surfaces.  The  forward 
transform  formula  requires  a  bit  more  work  to  develop  then  simply  substituting  Yr  into  (7a) 
for  Ys  (as  was  done  for  the  inverse  transform)  as  this  would  not  satisfy  the  consistency 

requirement   F[F_1(w)]  =  w.  To  derive  the  forward  transform,  both  sides  of  (12)  are 

alternately  multiplied  by  e±ipz ,  integrated  with  respect  to  z,  and  then  appropriate  linear 
combinations  are  taken.  The  result  of  this  process  yields  the  following  for  the  forward 
Fourier  transform: 

1 


u 


(x,p)=  ju(x,z)eipzdz  +  —Y^  \u{x^z)^  'pz^z  +  i \k\P^p{x^)^  ,  (13) 

0  *■  r\P)  0  0 


where  the  kernal  k[p,q)  in  the  last  integral  on  the  right  hand  side  is: 

rr(p)-rr{g) ,  i-rr{p)rr{qj 


k(p,q)  = 


+ 
q-p  q  +  p 


2*rr(P)  ' 


(14) 


For/?  =  q,  a  limit  can  be  performed  to  yield: 


k(p,p)  = 


-rip) 


+ 


i-r?(/?) 


2/7         2nrr{p) 


(15) 


where   r/(/?)    represents  the  derivative  of  the  rough  surface  reflection  coefficient.  In 
deriving  (13),  we  have  made  use  of  the  following  identity  from  Papoulis  [7]: 


0  P 


(16) 


where  8(-)  is  the  delta  function. 

Equation  (13)  may  be  rewritten  in  operator  form  by  denoting  the  last  integral  on  the 
right  hand  side  as,  Kcii  which  yields: 


u(x,p)  =  (l-iKc) 
=  ¥r[u{x,z)] 


oo  1  00 

u(x,z)eipzdz  +  — 7— Y    u\x,z)e~ipzdz 


(17) 


where  Fr  represents  the  forward  Fourier  transform  for  rough  surfaces  and  /  is  the  identity 
operator.  For  notational  convenience  we  define  the  quantity  P: 

P  =  {l-iKc).  (18) 
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III.  SOLUTION  PROCEDURE 

In  the  previous  chapter  we  presented  the  case  for  using  the  parabolic  equation  and 
split-step  algorithm  to  tackle  the  propagation  problem  of  vertically  polarized  waves  over 
rough  ocean  surfaces.  We  then  defined  a  consistent  forward  and  inverse  Fourier  transform 
pair  to  be  used  with  the  split-step  algorithm.  In  this  chapter  we  explain  the  solution 
procedure  using  this  transform  pair  to  determine  field  strengths  at  a  given  range. 

The  solution  procedure  begins  by  recalling  that  we  defined  our  transmitter  source  to 
be  an  omnidirectional  point  source  at  a  height  z  =  zx.  Numerically  this  initial  field,  w(x,z)  at 
range  x  =  0,  can  be  represented  by  a  delta  function.  To  determine  the  field  at  our  first  step 

downrange,  (x  =  nAx)  forn  =  1 ,  (17)  is  solved  for  u\x,p)  with  u(x,z)  equal  to  the  delta 
function.  For  this  initial  case,  (17)  becomes: 


u(x,p)  =  (l-iKc) 


e~'pz' 
eipz'+- 


(19) 


r,WJ  • 

Once  the  initial  u\x,p)  is  known  the  u  at  the  new  range  is  determined  by  propagating  in 


Ax 
->P~ 


free  space.  This  is  represented  by  the  term  e     2k°  in  (2).  Equation  (12)  can  then  be  solved 

for  the  field  u(x,z)  at  the  first  step  downrange.  A  second  phase  correction  is  then  applied  to 
this  field  to  compensate  for  variations  in  the  index  of  refraction.  This  phase  correction  is 
seen  as  the  e'*°(m-,)A,r  term  in  (2). While  this  phase  correction  is  suitable  for  most 
environmental  scenarios,  it  is  not  considered  adequate  for  rough  surfaces  when  severe 
refractive  index  gradients  are  present.  For  this  case  we  include  the  first  and  second 
derivatives  of  the  refractive  index  and  then  modify  (2)  as  shown  below.  The  first  and  second 
derivatives  were  determined  by  doing  a  cubic  spline  interpolation  of  the  given  refractivity 
index  data.  The  new  formula  that  includes  higher  derivatives  of  the  refractive  index  is: 
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■Z.       A     [(  ,\     ^'^l2 

*V"Axl  (.m-U 7 — 

u(x,z)  =  e 


[m"Ax]2 


-ip 


*b1 


?r[w(x0,z)] 


(20) 


where  Ax  =  x  -  x0 .  Once  a  solution  to  (20)  has  been  determined  we  can  now  take  another 
step  downrange,  (x  =  nAx)  for  n  =  2 ,  and  repeat  the  process  described  above.  This  range 
stepping  process  continues  until  the  field  is  determined  at  the  desired  distance  downrange. 

For  a  given  u(x,z)  (17)  is  solved  by  numerically  evaluating  the  quantity  in  the 
square  brackets  and  then  applying  the  inverse  operator  P"1  on  it.  The  numerical  solution  was 
obtained  using  MATLAB  code  to  compute  the  Fourier  transforms  by  means  of  an  N  point 
complex  FFT.  In  order  to  program  the  solution  to  (1 7)  the  continuous  operator  P  first  had  to 
be  put  into  a  discrete  form.  This  was  done  by  using  Simpson's  rule  with  weights,  Sn,  to 
discretize  the  last  integral  in  equation  (13)  to  yield: 


2 


Kcu[x,mAp)  «  Ap /  , SJu[x, nApikymAp, nAp) 


n=\ 


(21) 


=  Ku 


,  N     N) 
where  K  is  of  order  I  —  x  —  I  and  is  the  discrete  version  of  the  continuous  Kc.  The  upper 

N 
limit  in  the  summation  is  —  and  not  A^  because  the  integral  is  over  the  semi-infimte 

interval  [0,co)  ,  whereas  the  complex  FFTs  assume  limits  of  (-  oo,oo) .  The  lower  limit  is  1 
because  w(x,0)  =  0 .  The  elements  of  K  are: 


Kmn  =  ApSnk(mAp,nAp)  ,m,n  =  1,2,...  —  , 


(22) 


As  mentioned  above,  an  Appoint  FFT  is  used  for  the  computation  of  the  Fourier 
transform  pair.  Let  us  assume  that  the  various  quantities  are  band  limited  over  -/?max  <  p  < 
pmax,  and  that  the  transform  is  evaluated  atp  =  0,  A/?,  2Ap,  ...,  (N-\)  A/?,  where: 
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Ap  = 


2n 

NAz 


(23) 


(N     \ 
Positive  wavenumbers  occur  at  p  =  Ap,2Ap,...,\  —  -1  \Ap,  while  negative  wavenumbers 

(N     \       (N      \  ,         .  N 

occur  at  \—  +  \\Ap,\—  +  2\Ap,...,{N-\)Ap  .  The  value  —  Ap  corresponds  to  both  -pmax 

andpmax. 


Ap) 


Figure  3.  Wavenumber  Spectrum  p 

To  contain  the  computational  domain  vertically  in  the  physical  space,  up  to  a 
maximum  height  zmax,  and  to  bandlimit  the  signal  in  the  /?-space,  a  Harming  window  is 
employed  in  both  the  spatial  and  wavenumbers  domains  as  follows: 
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h(n)  = 


3N 
0<n<T 

Ann      3N  N 


sin 


<n<  — 

N  '      8  2 


Note  that  the  Harming  window  forces  a  gradual  rolloff  to  zero  of  the  data  set  as  can  be  seen 
in  Figure  4. 
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Figure  4.  Manning  Window 
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IV.  RESULTS 

In  Chapters  I  -  III  we  presented  theoretical  background  and  described  the  solution 
procedure  for  using  the  parabolic  equation  and  Fourier  split-step  algorithm  to  solve  the 
problem  of  propagation  of  vertically  polarized  waves  over  rough  ocean  surfaces.  In  this 
chapter  we  present  numerical  results  for  a  variety  of  propagation  scenarios.  We  will 
investigate  the  effects  on  propagation  by  varying  parameters  such  as  wind  speed,  step  size, 
frequency,  height  of  the  upper  boundary  and  refractivity  index  profile.  Throughout  this 
chapter  we  present  plots  showing  Propagation  Factor  (PF),  in  dB,  versus  receiver  height  or 
distance  downrange,  where  PF  is  defined  as  excess  signal  over  free  space.  For  a  point 
source,  the  propagation  factor  can  be  computed  from: 

PF  =  101og(H2x^0)  ,  (24) 

where  jc  is  the  downrange  or  horizontal  distance  from  the  transmitter,  u  is  the  field,  and  k0  is 
the  operating  wavelength  of  the  propagating  signal.  A  positive  (negative)  value  of 
propagation  factor  indicates  a  gain  (loss)  with  respect  to  propagation  in  free  space. 

A.    CASE  1  -  STANDARD  ATMOSPHERE 

We  begin  our  numerical  results  by  running  a  test  case  against  known  data  to  validate 
our  PE  formulation.  We  look  at  propagation  in  standard  atmosphere  where  the  refractivity 

M  =  ( m  - 1)  x  1 06  is  given  by: 

M(z)  =  (340  +  0.118z)  ,  (25) 

where  z  is  the  height  in  meters.  The  transmitter  is  at  a  height  of  z,  =  30m,  horizontal  step 
size  is  Ax  =  200m  and  the  ground  constants  are  given  by  sr  =  80  and  a  =  4  Seimens  per 
meter  (S/m).  Figure  5  shows  the  refractivity  profile  for  the  case  of  standard  atmosphere. 
Figure  6  shows  Propagation  Factor  versus  Receiver  Height  at  a  range  of  40  km  with  an 
upper  boundary  Zmax  =  512m,  frequency  of  operation  =  3  GHz  and  N=  512.  This  lobing 
pattern  compares  very  favorably  with  that  given  by  Kuttler  and  Dockery  [3].  Figures  (7-10) 
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illustrate  the  effects  of  increasing  wind  speed  and  hence  increasing  surface  roughness,  for 
various  heights  of  the  upper  boundary.  Note  that  while  the  value  of  the  upper  boundary 
changes  between  Figures  7,  8  (150m)  and  Figures  9,  10  (300m),  the  vertical  step  size  Az 
remains  the  same  (0.146m).  This  is  so  because  when  the  height  of  the  upper  boundary  was 
doubled  the  number  of  points,  N,  in  our  FFT  also  doubled.  For  these  figures,  frequency  = 
10  GHz,  transmitter  and  receiver  are  both  at  a  height  of  25m,  and  Az  =  0.146m.  From 
Figures  7  and  8  one  can  see  that  as  wind  speed  increases  the  excursions  of  the  PF  in  the 
interference  region  is  reduced  and  at  longer  ranges  the  PF  curve  no  longer  decays  smoothly. 
The  departure  of  propagation  factor  from  monotonous  decay  for  large  ranges  is  due  to 
numerical  reflections  from  the  upper  boundary.  Figures  9  and  10  have  an  upper  boundary  of 
300m.  We  can  see  by  increasing  Zmax  (and  keeping  the  step  size  Az  constant)  the  numerical 
reflections  from  the  upper  boundary  do  not  affect  the  decay  of  the  PF  until  further 
downrange. 

B.         CASE  2  -  TRI-LINEAR  DUCT 

Our  second  case  involves  examining  propagation  where  the  index  of  refraction  is 
characterized  by  a  tri-linear  duct  as  shown  in  Figure  1 1 .  Numerically,  refractivity  is  given 
by: 


M(z)  = 


340  +  0.118z         0<z<135, 
499.03-l.06z       135<z<150: 
322.33  +  0.118z    z>150, 


where  z  is  in  meters.  In  this  example  the  frequency  of  operation  is  3  GHz,  wind  speed  Ws  = 
0,  transmitter  height  =  30m,  and  Zmax  =  512m.  We  choose  N  =  512,  resulting  in  a  vertical 
step  size  Az  =  2m,  and  choose  a  horizontal  step  size  Ax  =  200m.  Figure  12  shows  the 
Propagation  Factor  versus  Receiver  Height,  at  a  range  of  40  km.  As  for  the  case  of  a 
standard  atmosphere,  the  propagation  characteristics  for  the  tri-linear  duct  are  in  agreement 
with  Kuttler  and  Dockery  [3]. 
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C.         CASE  3  -  EVAPORATION  DUCT 

The  next  refractivity  profile  we  investigate  is  that  of  an  evaporation  duct  as  depicted 
in  Figure  13.  The  numerical  values  for  the  evaporation  duct  are  shown  in  Table  1 . 
Table  1.  Refractivity  Data  For  Evaporation  Duct 


Height  (ra) 

M(z) 

0.000 

340.00 

0.135 

323.00 

0.223 

321.76 

0.368 

320.53 

0.607 

319.31 

1.000 

318.11 

1.649 

316.94 

2.718 

315.83 

4.482 

314.80 

7.389 

313.91 

12.182 

313.26 

20.000 

312.99 

33.115 

313.38 

54.598 

314.81 

90.017 

317.99 

148.413 

324.04 

165.000 

325.76 

300.000 

339.745 

Figures  14  and  15  illustrate  that  by  wisely  choosing  the  value  of  the  upper 
boundary,  accurate  data  can  be  obtained  where  it  may  not  have  been  otherwise  possible. 
Figure  14  shows  propagation  in  an  evaporation  duct  with  Zraax=75m,  N=512,  transmitter  and 
receiver  height  =  25m,  and  frequency=10  GHz.  At  the  longer  ranges  the  effects  of  the  upper 
boundary  and  the  evaporation  duct  prevent  the  PF  from  decaying  in  a  stable  manner.  By 
increasing  the  upper  boundary  to  150m  and  N=1024  and  keeping  all  others  parameters  the 
same,  Figure  1 5  shows  a  stable  PF  can  be  achieved  even  at  the  longer  ranges. 


D. 


EVAPORATION  DUCT  WITH  ROUGH  SEA  SURFACE 


In  this  example  we  consider  the  same  evaporation  duct  as  above,  but  now  we  add 
the  effects  of  a  rough  ocean  surface.  The  wind  speed  considered  is  10  m/s  which  results  in  a 
wave  height  of  0.51m  from  (11).  Figure  16  shows  both  the  propagation  factors  for  a  signal 
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in  the  evaporation  duct  with  no  wind  Ws  =  0  m/s  and  for  Ws  =  10  m/s.  Other  data  are  Zmax  = 
150m,  frequency  =  10  GHz  and  N  =  1024.  The  effects  of  a  rough  ocean  surface  are  more 
pronounced  than  the  case  of  standard  atmosphere.  Specifically  we  see,  (i)  a  reduction  of  the 
specular  component  as  seen  by  the  smaller  excursions  of  the  rough  sea  PF  in  the 
interference  region,  and  (ii)  increased  losses  for  the  rough  sea  PF  over  the  smooth  sea  PF  at 
greater  distances.  Figure  17  gives  the  results  for  PF  vs.  Receiver  Height  for  the  evaporation 
duct  at  a  range  of  100  km  for  wind  speeds  of  0  m/s,  10  m/s  and  20  m/s.  The  upper 
boundary  is  Zmax  =  300m,  frequency  of  operation  =  10  GHz,  N=2048,  and  transmitter  height 
z,  =  25m.  This  figure  illustrates  the  overall  reduction  in  PF  as  surface  roughness  increases, 
however  it  also  shows  significant  variation  in  the  value  of  PF  with  respect  to  height  due  to 
the  effect  of  the  evaporation  duct. 

E.         CASE  4  -  SURFACE  DUCT 

The  last  example  to  be  considered  is  the  case  of  propagation  in  a  surface  duct.  This 
refractivity  is  depicted  in  Figure  1 8  and  its  numerical  form  is  given  by  : 

|350-0.335z  0<z<45.7 

M^z>=  [329.36  +  0.1 164z    z>45.7 

Figure  19  shows  the  Propagation  Factor  vs.  Range  out  to  a  range  of  100  km  for  a  signal 
traveling  in  the  surface  duct  over  smooth  and  rough  ocean  surfaces  (Ws  =  0  m/s  and  Ws  = 
10  m/s).  The  frequency  of  operation  is  10  GHz,  transmitter  and  receiver  height 
z,  =  25m  =  zr ,  and  an  upper  boundary  Zmax  =  150m.  As  in  the  previous  example,  Figure  19 
illustrates  that  the  rough  ocean  surface  generally  decreases  the  propagation  factor  as  wind 
speed  increases.  The  last  example  considered  is  Propagation  Factor  vs  Receiver  Height  at  a 
range  of  100km  for  the  surface  duct  with  an  upper  boundary  Zmax  =  300m.  The  results  are 
shown  in  Figure  20  with  all  parameters  remaining  the  same  as  the  previous  case  except  here 
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N  =  2048.  Once  again  surface  roughness  has  reduced  the  propagation  factor,  and  the  effect 
of  the  surface  duct  has  resulted  in  variations  of  the  PF  as  receiver  height  is  changed. 
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Figure  5.  Refractivity  for  Standard  Atmosphere. 
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Figure  6.    PF  vs.  Receiver  Height  for  standard  atmosphere  at  range  of  40  km. 

Transmitter  height  zt  =  30m,  wind  speed  Ws  =  0,  frequency  =  3  GHz,  N 
512,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  7.    Smooth  and  rough  sea  (Ws  =  10  m/s)  PF  vs.  Range  for  standard 

atmosphere  out  to  a  range  of  100  km  with  upper  boundary  Zmax  =  150m. 
Transmitter  height  zt  =  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz, 
N  =  1024,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  8.    Smooth  and  rough  sea  (Ws  =  20  m/s)  PF  vs.  Range  for  standard 

atmosphere  out  to  a  range  of  100  km  with  upper  boundary  Zmax  =  150m. 
Transmitter  height  zt  =  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz, 
N  =  1024,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  9.    Smooth  and  rough  sea  (Ws  =  10  m/s)  PF  vs.  Range  for  standard 

atmosphere  out  to  a  range  of  100  km  with  upper  boundary  Zmax  =  300m. 
Transmitter  height  zt  =  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz, 
N  =  2048,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  10.  Smooth  and  rough  sea  (Ws  =  20  m/s)  PF  vs.  Range  for  standard 

atmosphere  out  to  a  range  of  100  km  with  upper  boundary  Zmax  =  300m. 
Transmitter  height  zt  =  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz, 
N  =  2048,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  11.  Refractivity  for  Tri-Linear  Duct 
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Figure  12.  PF  vs.  Receiver  Height  for  tri-linear  duct  at  range  of  40  km.  Transmitter 
height  zt  =  30m,  wind  speed  Ws  =  0,  frequency  =  3  GHz,  N  =  512,  vertical 
polarization  and  omnidirectional  antenna. 
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Figure  13.  Refractivity  for  Evaporation  Duct. 
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Figure  14.  PF  vs.  Range  for  evaporation  duct  to  a  range  of  100  km  with  upper 

boundary  Zmax  =  75m.  Transmitter  height  z,  =  25m,  receiver  height  zr  = 
25m,  frequency  =  10  GHz,  wind  speed  Ws  =  0,  N  =  512,  vertical 
polarization  and  omnidirectional  antenna. 
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Figure  15.  PF  vs.  Range  for  evaporation  duct  to  a  range  of  100  km  with  upper 

boundary  Zmax  =  150m.  Transmitter  height  z,  =  25m,  receiver  height  zr  = 
25m,  frequency  =  10  GHz,  wind  speed  Ws  =  0,  N  =  1024,  vertical 
polarization  and  omnidirectional  antenna. 
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Figure  16.  Smooth  and  rough  sea  (WS  =  10  m/s)  PF  vs.  Range  for  evaporation  duct  to 
a  range  of  100  Ian  with  upper  boundary  Zmax  =  150m.  Transmitter  height  zt 
=  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz,  wind  speed  =  0  m/s, 
N  =  1024,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  17.  PF  vs.  Range  for  evaporation  duct  to  a  range  of  100  km  with  wind  speed  = 
0  m/s,  10  m/s,  20  m/s.  Upper  boundary  Zmax  =  300m,  transmitter  height  zt 
25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz,  N  =  2048,  vertical 
polarization  and  omnidirectional  antenna. 
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Figure  18.  Refractivity  for  Surface  Duct. 
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Figure  19.  Smooth  and  rough  sea  (WS  =  10  m/s)  PF  vs.  Range  for  surface  duct  out  to  a 
range  of  100  km  with  upper  boundary  Zmax  =  150m.  Transmitter  height  zt  = 
25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz,  wind  speed  Ws  =  0 
m/s,  N  =  1024,  vertical  polarization  and  omnidirectional  antenna. 
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Figure  20.  PF  vs.  Receiver  Height  for  evaporation  duct  at  a  range  of  100  km  with 

wind  speeds  Ws  =  0  m/s,  10  m/s  and  20  m/s.  Upper  boundary  Zmax  =  300m. 
Transmitter  height  zt  =  25m,  receiver  height  zr  =  25m,  frequency  =  10  GHz, 
N  =  2048,  vertical  polarization  and  omnidirectional  antenna. 
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V.  CONCLUSIONS 

In  this  thesis,  a  method  of  predicting  radiowave  propagation  of  vertically  polarized 
waves  over  rough  ocean  surfaces  was  implemented  and  tested.  An  efficient  solution  was 
made  possible  by  only  considering  one  way  propagation,  thus  allowing  us  to  use  the  split- 
step  parabolic  equation  method  to  tackle  our  problem  instead  of  the  much  more 
computationally  intensive  Helmholtz  equation.  Our  method  incorporates  the  contributions 
from  direct  and  reflected  waves  at  any  distance  downrange  from  the  transmitter,  but  ignores 
the  effects  of  backscattering.  This  allowed  us  to  develop  a  range  stepping  technique  to 
determine  the  strength  of  the  propagating  waves.  The  fact  that  the  split-step  technique 
involves  Fourier  kernels  (plane  waves)  made  it  possible  to  incorporate  the  effects  of  sea 
surface  roughness  directly  into  the  spectral  domain  through  the  use  of  a  rough  surface 
reduction  factor. 

The  main  purpose  of  this  thesis  was  to  modify  the  equations  developed  by 
Janaswamy  [1]  so  the  propagation  of  vertically  polarized  waves  could  also  be  solved  by  the 
split-step  PE  method.  Numerical  results  for  the  case  of  vertical  polarization  show  great 
similarity  with  those  described  in  that  report  for  the  various  refractivity  profiles  examined. 
The  effects  of  the  upper  boundary  layer  are  a  significant  factor  in  the  accuracy  of  our  model 
and  must  be  taken  into  consideration  when  interpreting  data.  The  user  must  take  into 
consideration  the  value  of  the  upper  boundary,  surface  roughness  and  vertical  step  size  as 
determined  by  the  number  of  points,  N,  used  in  the  FFT,  to  ensure  accurate  results  are 
obtained  with  this  model. 

The  model  presented  in  this  thesis  is  applicable  to  propagation  predictions  for  areas 
such  as  communications  and  radar.  It  is  a  useful  tool  for  those  designing,  developing,  and 
deploying  new  systems  and  for  those  that  may  be  analyzing  systems  already  in  use. 
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